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MICROASPERITY MODEL FOR ELASTOHYDRODYNAMIC LUBRICATION 
OF A SPINNING BALL ON A FLAT SURFACE 
by Charles W. Allen* and Erwin V. Zaretsky 
Lewis Research Center 

SUMMARY 

A microasperity elastohydrodynamic (EHD) model was developed for a ball spinning 
without rolling on a flat surface. The theory model was developed by applying to a sinu- 
soidal surface model, existing dry- contact deformation theory and point- contact elasto- 
hydrodynamic theory. The computed torques based on this model were compared with 
experimental torques using the NASA spinning torque apparatus. 

There is fair agreement between the theoretical and experimental results below 
7 2 

300 000-psi (207x10 -N/m ) maximum Hertz stress although the exponential form of the 
torque-stress curves differ. Theory would indicate that microasperity EHD lubrication 
would be most effective where the nominal Hertz stress is low and the sliding velocity is 
high. Both of these conditions are satisfied at the outer region of the contact circle. 

This, however, is the region where the first asperity contact does occur if a test is run 
for long enough duration. Based upon the model presented, microasperity EHD lubrica- 
tion does not appear to be a primary force in maintaining separation between the opposing 
surfaces. 


INTRODUCTION 

There are several factors that can affect ball spinning torque and rolling friction in 
an angular- contact ball bearing. Some of the parameters which have been reported to 
affect ball-spinning torque are ball-race conformity, contact stress, and lubricant type 
(refs. 1 to 5). Other factors that may affect the spinning torque are surface finish, spin- 
ning speed, lubricant viscosity, and other lubricant properties. 

Much work has been performed to define the lubricant and its behavior in a noncon- 
forming groove where spinning without rolling occurs (refs. 5 and 6). In reference 6 
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an elastohydrodynamic (EHD) model was developed. However, it was assumed that some 
type of lubricant film existed within the center region of contact (fig. 1) within which it 
was not possible to generate an EHD film. There was good correlation between the anal- 
5 rtical and experimental results. However, there were several factors not considered in 
the analysis which may or may not be of secondary importance. One of these factors is 
the effect of the surface roughness which may introduce a microasperity lubrication 
effect. Microasperity lubrication, whether it be hydrodynamic or elastohydrodynamic, 
in which the deformation and pressure-viscosity effects must be considered, is a phen- 
omenon which, in the general case, must include the shape, size, deformation and dis- 
tribution of the asperities at the lubricant interface and the interaction between individual 
asperities. Very little work on the theory of microasperity elastohydrodynamic lubrica- 
tion has been reported. 


Inscribed circle defined 



(a) Contact region for ball in nonconforming groove. 



(b) Contact region for ball on flat surface. 

Figure 1. ~ Contact geometries for ball on nonconforming surface. 


Some work has been reported (refs. 7 and 8) in microasperity lubrication in which 
the viscosity is assumed to be constant and deformation is neglected. Both of the refer- 
enced works consider the asperities to be cylinders distributed uniformly on a plane and 
in nominal contact with another smooth plane. The lubrication of sinusoidal asperities 
by a viscoelastic fluid was considered in reference 9. Elastohydrodynamic theory has 
been applied to the collision of two individual asperities (ref. 10). The asperities are 
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considered to be of parabolic shape and long; therefore, the problem is reduced to a 
two-dimensional analysis. 

The objectives of the work reported herein were (1) to develop a model to describe 
microasperity elastohydrodynamic lubrication theory for a ball spinning without rolling 
on a flat surface and (2) to compare the analysis with experimental results. The method 
of analysis considers first an asperity surface model for dry-contact. Subsequently, an 
elastohydrodynamic solution is undertaken using point-contact theory for the microas- 
perities. The dry- contact and EHD solutions were then combined into a microasperity 
elastohydrodynamic theory as suggested in reference 11. The results of this theory were 
compared to the experimental results reported in reference 2. 


SYMBOLS 

2 2 

A area, in. (m ) 

a radius of contact circle, in. (m) 

d distance between two surfaces, in. (m) 

E modulus of elasticity, psi (N/m ) 

2 

H Vickers hardness, psi (N/m ) 

K thermal conductivity, Btu/(sec)(ft)(°F) ; (W/(m)(‘^C)) 

I average distance between asperities, in. (m) 

M spinning torque , lb-in. (N-m) 

o 

2 

p pressure, psi (N/m ) 

2 

Pj^ pressure at which viscosity exponent changes, psi (N/m ) 

p^ average pressure over typical asperity, psi (N/m ) 

Q heat generation, Btu/sec (W) 

q heat flux, Btu/(ft^)(sec) (W/m^) 

R radius of ball, in. (m) 

r radius of elemental annulus , in. (m) 

T temperative, °F (°C) 

t time, sec 

W normal load, lb (N) 

V entraining velocity (vector sum of surface velocities), in. /sec (m/sec) 
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Ill 


2 

a,l3 pressure- viscosity exponents, psi (N/na ) 

2 2 

K dtffusivity, in /sec (m /sec) 

2 2 

jU absolute viscosity, reyns, Ib-sec/in. (N-sec/m ) 

2 2 

absolute viscosity at entry to contact zone, reyns, Ib-sec/in. (N-sec/m ) 
p radius of curvature, in. (m) 

a composite surface roughness (rms), ju in. (p cm) 

CTj^,ct 2 surface roughness (rms) for surfaces 1 and 2, /i in. (p cm) 
ct' amplitude of sinusoidal surface model, p in. (p cm) 

T shear stress, psi (N/m ) 

u) angular velocity, rad/sec 


SURFACE ROUGHNESS 

The surface roughness models considered herein were developed for the case of dry- 
contact and normal approach with no sliding. Sliding under unlubricated conditions would 
cause shearing and plastic deformation of the asperities. However, if a thin film can be 
shown to exist between asperities, the shear will take place within this film and the nor- 
mal approach model would be applicable for a first approximation (ref. 12). 

Some of the proposed asperity models are shown in figure 2. The models protrayed 
may be treated as two-dimensional in which the hemispheres become cylinders, the cones 
become wedges, etc. They also may be treated three dimensionally. 

In reference 13, a sophisticated model is considered. This model consists of hemi- 
spherical asperities which have superimposed on them other much smaller hemispherical 
asperities. The distribution of asperity heights is usually close to Gaussian as reported 
in reference 14. 

Another very important aspect of the contact of rough surfaces is the mode of defor- 
mation. Frequently, the deformation is assumed to be elastic or ideal plastic although 
plastic hardening or elastic-plastic are sometimes considered, as reported in refer- 
ence 15. It is also shown in reference 15 that, if an exponential distribution is used in- 
stead of the more realistic Gaussian distribution, the true area of contact is directly pro- 
portional to the load. This is not true for the Gaussian distribution and a particular 
mode of deformation must be assumed. A true Gaussian distribution requires a finite 
probability of some plastic and elastic deformation existing together. After repeated 
contacts have been made, however, the high asperities have all yielded. During subse- 
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(a) Cylinders. 



(b) Hemispheres. 



(c) Cones. 



(d) Cones with spherical caps. 



(e) Ellipsoids. 



(f) Sinusoid. 

Figure 2. - Surface models. 


quent loading, the deformation will be essentially elastic with a load close to that re- 
quired for plastic deformation. 

In all surface contact work it is necessary to distinguish between the nominal and 
real areas of contact. The following definitions from reference 16 will be used through- 
out this report: 

(a) Nominal area of contact - the apparent area of overlap of the contacting solids 

(b) Real area of contact - the sum of the separate microscopic areas at which the 

microasperities are in physical contact 

The distribution of load over the nominal contact area will differ from the Hertzian 
distribution but in reference 17 it is shown that, at high loads, the overall pressure dis- 
tribution approaches the Hertzian shape. 

Much of the previous work on the contact of rough surfaces has considered the con- 
tact between a rough surface and a smooth surface. In reference 15 it is shown that, 
from a statistical viewpoint, the contact of two rough surfaces is the same as that of a 
rough surface and a smooth surface, provided the roughness on the single rough surface 
is taken as the composite value given by: 
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a = y aj + ( 7 ^ 


( 1 ) 


where: 

a composite roughness, (rms) 

roughness of surface 1, (rms) 

Ug roughness of surface 2, (rms) 

For microasperity elastohydrodynamic lubrication, it is the relative distance, d, 
between the two surfaces which is important (fig. 3). Therefore, the single rough sur- 
face with a composite roughness as given above, is deemed to be a sufficiently accurate 
representation of the real situation. 




Figure 3. - Surfaces in contact. 
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The surface model which is the most representative of typical bearing surfaces is 
the sinusoidal model and this will be used in the analysis. A typical sinusoidal asperity 
sliding against a smooth plane is shown in figure 4. In accordance with the relation be- 
tween the amplitude and rms values for a sinusoidal waveform, the ’’amplitude” of this 
sinusoidal form is taken as yT2o and the ’’wave length” i as the average distance be- 
tween asperities. The radius of curvature p at the peak is then given by 



4n^a’ 


( 2 ) 


where 

I average distance between asperities 

and 


O’’ = V2cr 

The load carried by an individual asperity will be computed on the basis that it is 
close to that at which plastic deformation occurs (ref. 14). In reference 18, curves are 
given which indicate that for a spherical contact, the onset of plastic flow occurs when 
the average Hertz pressure reaches about 40 percent of the Vickers Hardness Number. 
This will be taken as the average Hertz pressure over a typical asperity and is given by: 

p^ = 0.4H (3) 


average pressure over typical asperity 

Vickers hardness converted to same units as the average (Hertz) pressure 

If the nominal contact between the ball and flat plate is now considered, the true 
area of contact at a radius r from the center of the nominal circle of contact, radius a, 
may be determined by considering an elemental annular ring of width dr (fig. 5) and 
realizing that the nominal force, considering a Hertzian distribution, and the actual force 
applied through the asperities, are equal. 

p^ dA = p(2jrr dr) 


where 

Ph 

H 


where 
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dA 


elemental area of contact over aseprity peaks at radius r, therefore 


rfA ^ 2?rrp dr 

Pa 


But for a Hertzian contact: 


P = 


1. 5 W 


Tta, 



( 4 ) 



Figure 5. - Elemental annulus showing areas over which 
force Is applied through asperities. 


The radius of the nominal contact region a is given in reference 19 as 


a - 1.109 



( 5 ) 


where 


W 

normal load 

R 

ball radius 

E 

modulus of elasticity 

a 

radius of nominal contact region 

and 
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( 6 ) 


dA = 



dr 



ELASTOHYDRODYNAMIC LUBRICATION AT POINT CONTACTS 

Considerable experimental and some analsdical work has been performed on point- 
contact lubrication (refs. 20 to 23). An approximate expression for the film thickness 
under point- contact conditions is given in reference 21 as 

h = 0.67(a]U^jV)^/^p^/^ (7) 


where 

h film thickness 

a pressure -viscosity exponent 

viscosity of lubricant at entry to contact zone 
V entraining velocity (vector sum of surface velocities) 

p radius of curvature of point contact 

While there exists more sophisticated formulas to calculate film thickness, the 
above simple formula is deemed adequate as a first approximation for a microasperity 
elastohydrodynamic theory. 


SPINNING TORQUE 

For an individual asperity as shown in figure 6, the film thickness is given by equa- 
tion (7). Now, considering all the asperities contained within the annulus of radius r 
and width dr, the total area of contact is given by equation (6). If the fluid is assumed 
to behave in a Newtonian manner and the velocity gradient is assumed to be linear, then 
the shear stress t in the fluid on the top of the asperities is given by: 

^ wr 

T = p — 

h 
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where 

H viscosity within lubricating film 

0 ) angular velocity of ball spin about axis perpendicular to contact surface 


Asperity lubricating 



Figure 6. - IVtode of asperity lubrication. 


The elemental torque imposed by the microasperity film within this annulus is then 

dM„ - rr dA 

O ' 

( 8 ) 

h 


where 


M = Torque 
s 

In the above equation ju is constant because the pressure over the asperity peaks is as- 
sumed constant and, according to the analysis presented in appendix A, isothermal con- 
ditions may be assumed to prevail for tests of 15 to 30 seconds duration. For purposes 
of the model the effect of the fluid lying in the valleys between the asperities is neglected. 

The total spinning torque is thus the integral of equation (8) over the whole region of 
contact. 


M 


s 


= O) 



dA 


Substituting expressions for dA and h from equations (6) and (7), respectively, the 
spinning moment is obtained as: 
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M 


s 




I ^2/3 l/3„2„ 
(aii^) ' p ' a 



dr 


Integrating this numerically and substituting the value of a from equation (5) the spinning 
torque is obtained as 


0. 89 


RESULTS AND DISCUSSION 


Values of spinning torque were calculated using equation (9) for a ball spinning with- 
out rolling on a flat surface using a synthetic paraffinic oil without additives as the lubri- 
cant for the following conditions : 


Speed, rpm (rad/sec) 1050 (110) 

Load, lb (N) 5 to 60 (22 to 267) 

Maximum Hertz stress, psi (N/m^) 160 000 to 370 000 (llOxio"^ to 260x10"^) 

Surface roughness of ball, p in. (p cm) rms 2 (5. 1) 

Surface roughness of groove, p in. (p cm) rms 3 (7.6) 

Typical number of asperities per in. (per cm), ref. 16 1000 (394) 

Hardness of ball and flat, Rc 63 


Using the above values the composite roughness from equation (1) is found to be 

approximately 4 microinches (10.2 microcentimeters) rms, and the average asperity tip 

radius is approximately 0.005- inch (0.012 cm). From equation (3) the average pressure 

7 2 

over a typical asperity is 440 000 psi (306x10 N/m ). 

The pressure-viscosity relation for a synthetic paraffinic oil (ref. 6) is shown in 
figure 7 and given as : 


a Pa 

= P^e p^ ^ Pi 

(10a) 


(10b) 
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II III II II 


where 

o! pressure-viscosity coefficient; 9.2x10”^ (psi)"^, (1. 33x10"^ (N/m^)"^) 

/3 pressure- viscosity coefficient; 5x10“® (psi)“^, (0.72x10^ (N/m^)"^ 

Pj^ critical pressure; 55 000 psi (37x10*^ N/m^) 

ji viscosity at entry to asperity contact; 6x10”® reyns (0.4 N-sec/m^) 

Values of torque calculated from equation (9) are shown in figure 8 and compared to 
experimental results from reference 2. 


lU 


4000x10:5 


-P“ 5x10"^ (psi)"l [0.72x10’’ (N/m^l’^] 


-a = 9. 2x10’’ (psi)’l [l. 33x10’^ (N/m)’l] 


0 40 80 120 160 200 240 280 

Maximum Hertz stress, ksi 


0 20 40 60 80 100 120 140 160 180x10^ 

Maximum Hertz stress, N/m^ 

Figure 7. -Theoretical pressure-viscosity relation for syn- 
thetic paraffinic oil at 83° F (302 K) (from ref. 6), 


Theoretical/ ^ 

/ / 

/ / 

/ /Experimental (ref, 2) 


150 200 250 

Maximum Hertz stress, ksi 


Maximum Hertz stress, N/m"^ 

Figure 8. - Torque computed from microasperity EHD analysis compared with that experimentally 
obtained. 
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A comparison of the theoretical and experimental results at stress levels below 
300 000 psi (207x10 N/m'^) maximum Hertz stress shows them to be in agreement al- 
though the exponential forms of the curves differ. The effect of the fluid filling the 
valleys was ignored in the analysis. However, the torque contribution from this fluid 
may not be negligible. The exact contribution is dependent upon the viscosity of the fluid 
within the valleys. This in turn depends upon the pressure, which, for this analysis, 
was assumed to be atmospheric. The pressure could, however, be much greater. Any 
increase in the viscosity within the valleys will increase the film thickness over the 
peaks since the latter is dependent upon the viscosity of the fluid entering the real contact 
zone. This increased film thickness would reduce the shear stress over the asperities. 
Thus, the effect of the shear within the valleys may be decreased by the reduced shear 
force over the tips of the asperities. 

The area of real contact over the tips of the microasperities is given as a percentage 
of the nominal contact area in figure 9. The derivation of this relation is presented in 
appendix B. From figure 9 it may be seen that at stresses greater than 300 000 psi 
(207x10 ' N/m"^) over 65 percent of the nominal contact area is in actual microasperity 
contact near the center of the nominal contact circle. This traps fluid in "lakes” be- 
tween the asperities. The resulting squeeze effect would tend to increase the separation 
of the surfaces. Any increase in film thickness would tend to reduce the torque below the 
predicted value and this would be more evident at higher loads. Thus, the net effect 
would be to flatten the torque-stress curve. 



Figure 9. - Asperity contact area as percent of 
nominal area. 
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The film thickness between the tips of the asperities and the opposing surface as cal- 
culated by equation (7) is about 0.2 microinch (0. 5 Mem) (50 A). This is only a few 
molecules thick. Therefore, the validity of this analysis depends upon the bulk viscosity 
of the lubricant remaining unchanged for films a few molecules thick. Direct experi- 
mental verification of the viscosity of such thin films appears to be impossible at the 
present time. 

Based upon the assumptions made herein, microasperity EHD lubrication would be 
most effective where the nominal Hertz stress is low and the sliding velocity is high. 

Both of these conditions are satisfied at the outer region of the contact circle (fig. 10). 
This, however, is the region where the first metal to metal contact occurs during test. 
Therefore, microasperity EHD lubrication does not appear to be a primary force in main- 
taining separation between the opposing surfaces. 



0 . 005 . 010 . 015 . 020 

Distance from center of nominal contact circle, cm 

Figure 10. - Theoretical film thickness over microasperities as function of distance from center of nominal 
contact circle (eq. (7)). 


SUMMARY OF RESULTS 

A microasperity elastohydrodynamic model was developed for a ball spinning without 
rolling on a flat surface. The analytical results were compared with experimental re- 
sults previously obtained. The model assumed sinusoidal asperities model and com- 
bined the effects of dry contact deformation, point- contact deformation, and point- contact 
elastohydrodjmamic theory. The following results were obtained: 

1. There is agreement between the theoretical and experimental results at stress 
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n p 

levels below 300 000-psi (207x10 -N/m ) maximum Hertz although the exponential form 
of the torque-stress curves differ. 

2. Microasperity EHD lubrication would be most effective where the nominal Hertz 
stress is low and the sliding velocity is high. Both of these conditions are satisfied at 
the outer region of the contact circle. This, however, is the region where the first as- 
perity contact will occur during test. Based on the model presented herein, microas- 
perity EHD lubrication does not appear to be a primary force in maintaining separation 
between the opposing surfaces . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 30, 1970, 

126-15. 
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APPENDIX A 


ESTIMATE OF TEMPERATURE RISE AT BALL-PLATE INTERFACE 

As an example of the temperature rise at the interface consider the case of a maxi- 

7 2 

mum Hertz stress of 300 000 psi (207x10 N/m ) for which the experimental torque is 
0.018 pound- inch (0.00203 N-m). 

At a spinning speed of 110 radians per second the heat generated per unit time is 

Q = i = 2. 13X10'^ Btu/sec (0. 224 W) 

9336 

The area of the nominal contact circle is 

A = Tra^ = 1. 55x10"^ in. ^ (l.OxlO'”^ m^) 

If the heat flux is divided equally between ball and plate then the heat flux into either 
is 


q = Q = 0. 69 Btu/(in. ^)(sec) (1. 12x10® W/m^) 

A 

From reference 24 the heat flux into a semi-infinite solid over a circle of radius a, 
with no heat flux across any other boundary yields a maximum surface temperature rise 
T given by 


T = 



ierfc(O) - ierfc 



For steel: 


K = 5. 8x10"^ Btu/(sec)(in. )("^F) (4. 3 W/(m)(^C)) 

K= 1. 8x10"^ in. Vsec (1. 10x10"® m^sec) 

Using the foregoing equation and values, after 30 seconds the maximum temperature 
increase is 10*^ F (5.6^ C). The mean temperature increase would be less than this. 
Therefore the assumption of isothermal surfaces appears justified. 
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APPENDIX B 


NOMINAL AND REAL AREA OF CONTACT 


From equation (6) the elemental area of real contact at radius r is 


dA 



The nominal area of contact at this radius is 2nr dr. Therefore: 


Real contact 
Nominal contact 



This ratio expressed as a percentage is plotted in figure 9, 
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